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The Massive Soft Anomalous Dimension Matrix at Two Loops 
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We study two-loop anomalous dimension matrices in QCD and related gauge theories for products 
of Wilson lines coupled at a point. We verify by an analysis in Euclidean space that the contributions 
to these matrices from diagrams that link three massive Wilson lines do not vanish in general. We 
show, however, that for two-to-two processes the two-loop anomalous dimension matrix is diagonal 
in the same color- exchange basis as the one-loop matrix for arbitrary masses at absolute threshold 
and for scattering at ninety degrees in the center of mass. This result is important for applications 
of threshold resummation in heavy quark production. 



The infrared structure of perturbative amplitudes is 
relevant for a variety of hard-scattering processes, includ- 
ing the production of heavy particles, whether charged 
under QCD or not, and of high-p^ jets. Infrared en- 
hancements in these amplitudes are typically regularized 
by continuing to D dimensions, where they appear as 
poles in £ — 2 — D /2. Infrared poles, of course, do not 
appear in physical predictions for infrared safe quantities, 
in which they cancel after an appropriate sum over final 
states. Nevertheless, the all-order structure of infrared 
poles is of interest for exact fixed-order calculations 
for resummations of long-distance enhancements in cross 
sections and for supersymmetric gauge theories [1]. 

In this paper, we study the dependence of infrared 
enhancements at two loops on the masses of exter- 
nal partons. Our analysis is relevant to dimensionally- 
regularized amplitudes for 2 — > n processes (labelled by 

f : fi{Pi,ri) + h{V2,r2) 

/3(P3, fs) + • ■ ■ /n+2(Pn+2, r„+2) , (1) 

where the r,; represent color indices. At fixed, nonforward 
angles, amplitudes of this type factorize into functions 
representing short-distance (hard) dynamics, coUinear 
dynamics of external massless partons, and soft-gluon 
exchange between light and heavy partons. We can rep- 
resent such a factorized amplitude as [1, Q 



Hi A, 



Pi Q_ 



in the notation of Ref. [6], where the ket \A{p,e)) rep- 
resents X]l=i ^l{p,£) (cl)|-^.j, with {cl} some basis of 
color tensors linking the heavy and light parton lines at 
short distances. The jet functions for incoming and 
outgoing lightlike partons may be defined, for example, 
as the square root of the corresponding elastic form fac- 
tor 0. The soft matrix Sf, which is the subject of study 
here, describes the infrared behavior of color exchange 
between the external partons. It remains after the fac- 
torization of coUinear enhancements into the jet func- 
tions for those external partons that are massless, and 



in general mixes the color components found at short 
distance. We do not introduce jet functions for massive 
partons here. The resummation of coUinear logarithms 
associated with high-p^ massive partons can be treated 
by factorization as well [7], but this is not a goal of this 
paper. 

The soft matrix is determined entirely by a matrix of 
anomalous dimensions through 



Sf {Pi,as{fi),e) 
= P exp 



dp, _ ' 

M 



(3) 



where the velocities (3i are lightlike for light partons. 



0, and may be scaled to unit length for heavy 



particles, = 1. Integrals in the exponent are carried 
out in D > 4, using the D-dimensional running coupling, 
cts (/i, e). Equation ([3]), along with corresponding expres- 
sions for the jet functions, determines the infrared pole 
structure to all orders in perturbation theory for pro- 
cesses involving the wide-angle scattering of any number 
of massless and massive partons 0, H, [1] . 

The determination of the anomalous dimension ma- 
trix for an arbitrary process with massless and massive 
partons is equivalent to the renormalization of a set of 
color tensors that link the corresponding product of Wil- 
son lines at a point d, . Each Wilson line follows the 
velocity Pi of the corresponding parton, without recoil, 
from this point to infinity, either from the initial state 
or into the final state. These composite operators mix 
under renormalization in general, leading to the matrix 
structure shown in Eq. ([3]). The one- loop anomalous di- 
mensions for gluons and for both massless and massive 
quarks have been known for some time . At two 

loops, the matrices for any 2 —^ n process with massless 
lines satisfy the relation [ll| 



K 



(4) 



with T^gj the coefficient of (as/Tr)', and with K — 
Ca(67/18 - C(2)) - lOTpup/d. This is exactly the re- 
lation satisfied by the expansion of the cusp anomalous 
dimension [l2| , which generates the leading, double poles 
in the elastic form factor [l^, . 
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FIG. 1: Diagrams whose ultraviolet poles determine the soft 
anomalous dimension at two loops. The crossed vertex rep- 
resents the point at which the Wilson lines are linked. The 
straight lines represent eikonal propagators. Referring to the 
number of Wilson lines linked by gluons, in the text we refer 
to these as 3E diagrams (a-c) and 2E diagrams (d-f). 



Examples of the diagrams involved in the calculation 
of the two-loop anomalous dimension are shown in Fig. 
[TJ In momentum space, the propagators and vertices 
from Wilson lines are given by eikonal expressions Q. 
The corresponding two-loop corrections to the anomalous 
dimensions are found in the usual way from the two- 
loop UV single poles of these diagrams after one-loop 
renormalization. 

The result (|4]) for massless partons is a consequence of 
the vanishing of the single poles of those two-loop "3E" 
diagrams in which color is exchanged coherently between 
three eikonal lines in the figures. The arguments of Ref. 



ll| do not, however, generalize directly to massive Wil- 
son lines, with velocity vectors Pf ^ 0. While an ana- 

(2) 

lytic determination of V g would, of course, be desirable, 
numerical determination is also of interest, and is cer- 
tainly adequate to determine how far Eq. ([3]) generalizes 
to the production of massive particles. We provide the 
necessary analysis below, and show that when the (if are 
nonzero, Eq. ([¥]) no longer holds. A generalization of Eq. 
(HI), however, given by Eq. ([50)1 below, does holds for two- 
to-two processes for special momentum configurations. 

Much of our analysis will be carried out in position, 
rather than momentum space. In the following, we will 
take every parton as massive, and use the scale invariance 



of Wilson lines to set (if = \. Because we are calculating 
renormalization constants, we can carry out our analysis 
in Euclidean space. Indeed, a numerical result in Eu- 
clidean space is adequate to establish that the matrix 
does not follow Eq. (|4]) in Minkowski space. Otherwise, 
analytic continuation through Wick rotation would im- 
ply that the same result would hold in Euclidean space 
as well. 

We begin with the diagram. Fig. [T^, in which three 
eikonal lines are coupled by gluons that are linked at 
a three-gluon coupling [Tlj . In the configuration space 
evaluation of this diagram, we must integrate the posi- 
tion of the three-gluon vertex over all space. The three 
propagators each have one end fixed at this vertex and 
the other end fixed at a point Xiiii along the ith Wilson 
line. Each parameter Ai is integrated from the composite 
vertex to infinity. This diagram vanishes in Minkowski 
space for massless Wilson lines 

Suppressing color factors, we represent the 3E diagram 
Fig. [T^ as 



= d 



/'OO 



(5) 



Here /?/ = /32, /^s} denotes the set of three massive 
velocities of the lines to which the gluons attach, while 
the propagators and numerator factors of the integrand 
are given by a sum over six terms. 



(6) 



Each of these terms involves the derivative of one of the 
propagators, according to the usual gauge theory rules 
for the three-vector coupling, 

Vijk{x, Pi) = • Pj A(a; - XjPj) A(x - XuPk) 



X Pk ■ d^A{x - \iP,) 



(7) 



where A represents the position-space scalar propagator. 



A (a; - Xip^) 



r(i 



47r2- 



{x - XiPi) 



(8) 



We work in Feynman gauge. The contribution of this 
(scaleless) diagram to the anomalous dimension matrix 
is found from the residue of its simple ultraviolet pole. 
We note that all diagrams found from products of Wil- 
son lines are scaleless overall, and are defined by their 
renormalization constants 11 1 . 

At fixed for massive eikonals the A integrals in Eq. 
([5]) are all finite in four dimensions. After these integrals 
are carried out, the /3i-dependence enters only through 
the combination 



Pi ■ X 



(9) 
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and we can write 



(10) 



where Af{e) absorbs overall factors that are finite in the 
limit e = 0. To simplify our notation, in the following we 
normalize F^g so that Af{e) = 1. We recall that we have 
used the scale invariance of eikonal lines to set /3| = 1, 
and that / represents the set i,j,k. Each term ^ijk is 
now given by 



7jjfe(va;2,C/,£) = Pi ■ Pj f{x, f3j,e) f{x, (3k,e) 

X Pk ■ f{x,(3^,e) , (11) 

where the functions f{x, /3, e) are simply the integrals of 
the x-dependent factors of the propagators, 



dX 

(a;2 -2X(3-x + X^) 



l-e 



(12) 



After a change of variables to A' = A/va;^, the depen- 
dence on the variables and d factorizes, 



fix, Pi, e) 



l-2e 



dX' 



1 



2.., .0 (1-2A'0 + A'2) 

1 



— g(G,e)- 



(13) 



For the full expression, we also need the gradient of this 
function, which can be written as 



dUix,P,e) = 



1 



2-2e 



(2£-l)^5(C,e) 



x\ dg{C,e) 



(14) 



We note that this derivative is necessary to produce an 
overall x~'^ fall-off at infinity and a singularity at a; = 
0, corresponding to logarithmic infrared and ultraviolet 
behaviors. 

We next substitute the expressions for the functions 
/ in (|13p and their gradients into Eq. (jlip for the 
terms 7ijfe, to find 

ft • Pj 



(2e - l)Cfcff(0, e) + (-CfeC + ft ■ Pk) 



4-6e 



(15) 



In this expression, the first term in the square brackets 
is symmetric in the pair (i, j) and the third is symmetric 



in the pair [j, k) . The full nonvanishing contribution to 
Eq. (fTU]) is thus simply 



r ^ 

f^IHpi,^) ^ - d^x J2 'vkCkC 



ft • ft 



i j,/c— 1 



4-6£ 



X5(0,£)5(a,£)^^^. (16) 

Using the freedom to reintroduce dependence on the \f0l 
by demanding scale invariance, we can use this result in 
both Minkowski and Euclidean space to identify and iso- 
late the ultraviolet pole. It is now straightforward to 
show two important results that follow from the anti- 
symmetries built into Eq. p6|) . 

First, working in Minkowski space, we can readily con- 
firm the vanishing of F-^g for arbitrary massless ft. In 
this case, the function remains scale-invariant in the ft, 
although of course we cannot rescale by Pi = 0. Never- 
theless, the explicit form of g{(^, e) is 



1 1 

2^C 



(/32 = 0) 



(17) 



which, using (^(dg/dC,) — —g, immediately gives a vanish- 
ing integrand in Eq. (|16p by antisymmetry. It is inter- 
esting to note that, unlike the discussion in momentum 
space in Ref. [ll| , this proof of the vanishing of the three- 

gluon diagram, -Pjg , Fig. [T^, does not require a change 
of variables. 

In fact, the vanishing of Fig. [T^ extends to the case 
where only two of the three lines are massless [15]. Taking 
for definiteness = /3| = with /3| 7^ 0, and using Eq. 



U} in Eq. we find 

^4'^(ft'^) = - fd^'^—^ri^e 



X .9(C3,£) + C3- 



%(C3,£)^ /ft -ft ft -ft 



Ci 



C2 



(18) 



In this case, we follow [ll| and make a change of variables, 
using the light-like directions /3i,2 to define light-cone co- 
ordinates. To be specific, if we choose xi = Ci/ft 'ft £^nd 
X2 = C,2l Pi ■ Ps, we derive an integrand that is manifestly 
antisymmetric in xi and X2- We note that the momen- 
tum space method of [ll| also applies directly to the case 
when two of the three Wilson lines are massless, although 
this was not pointed out explicitly there. In both cases, 
the relevant change of variables exchanges two lightlike 
directions. This approach does not show that diagrams 
with a single massless line vanish identically, and indeed 
this seems unlikely, given that in this case there is only 
a single lightlike direction. 

Our second main result is that for both massive and 

f 21 

massless Wilson lines the function ^^JiPi ■ Pj) vanishes 
when any pair of the invariants are equal, say. Pi ■ P2 = 



4 



f3i ■ (33. We can show this by changing variables in Eu- 
chdean space from to r = Vx^ and the Q = (3i-x/Vx^. 
A straightforward calculation gives a form in which the 
overall scalelessness of the diagram is manifest in the ra- 
dial integral, 



dr 



dCidC2dC3 



X E ^^^■'^ ^« ■ 9iQ,e)g{Ck,e) • (19) 



Defining Wj = Pi ■ (3j, S,j = Q - WjCi, and Z3 = [(32 ■ P3 - 
W2Wz\/ \J\ — the explicit form of K{Qi) for arbitrary 
W2 and W3 is 



We now turn to the remaining 3E "double-exchange" 
diagrams, Fig. lb,c. These diagrams are given by four A 
integrations, two of which are along the Pj line. For ex- 
ample. Fig. lb may be written in the notation introduced 




H^-Cl){l-wl){l-zl^wl). 



- 

(20) 



The function K{C,i) is symmetric under the interchange 
of C2 and when W2 = W3. We easily check that when 
W2 = ^3 the remaining integrand in Eq. I|19p is antisym- 
metric under the exchange of C2 and Cs- The variables 
and are exchanged by a reflection in the two-sphere 
defined by any fixed value of Ci about the axis specified 
by the projection of P2 + Pa into this two-sphere. The 
integration region of (^2 and (.3 is therefore also symmet- 
ric. We conclude that the integral over ^2 and in Eq. 
(fTO)) vanishes when Pi ■ P2 = Pi ■ Pa- This relation holds 
as well in Minkowski space, which can be reached by an- 
alytic continuation at fixed ratios of the inner products 

Pi ■ P2/P1 -Ps. 

We do not have here an analytic form for the residue 

(2) 

of the UV pole of F^^' for generic values of the invari- 
ants Pi ■ Pj. A numerical analysis of Eq. (|16p . however, 
is particularly straightforward in Euclidean space. For 
this purpose, it is convenient to use D-dimensional polar 
coordinates, r, fto-i- The single overall ultraviolet pole 
in the scaleless integral (|16|) appears at r = 0, and the 
remaining three angular integrals at e = determine the 
residue of the pole in MS renormalization. These can 
be carried out readily, using the elementary form of the 
function g{C,0) in Eq. (fT^ . 



5(C,0) 



arccos ( 



(21) 



In Fig. [21 we show a plot of the residue of the UV pole 
of F^f for Pi ■ Ps = 0.5 in the Pi ■ P2/P2 ■ Ps-plane, 
suppressing the surface where the integral is negative for 
clarity. Notice the lines of zeros at Pi- P2 = /?2 • Ps and at 
P1-P2 — 0.5 and P2-P3 = 0.5. The peak towards P1-P2 
1 refiects an additional singularity in the integrand when 
two velocities are parallel. 



FIG. 2: Plot of the integral in ^ for Pi - Ps ^ 0.5 



above as 



Mb = CibJ^ie) {P,-P,) iPrPk) 



dX 



x/(A,,a/3j,A,£)/(A,,f,/3„/?fc,e), (22) 

where the functions / are defined as in Eq. (|13p with 
X = Xj,cPj and with c — a,b. 77{e) absorbs overall factors 
that are finite in the limit e = 0. In the following we 
set A/'(e) = 1. We have kept the overall color factor, 
represented by Cib- The variables C,i = Pi ■ x/Vx^ , are 
independent of the scale of x in Eq. ((22|) . so that all 
dependence on the \j^c variable is in the overall factor of 
x"^ = •^j,c> giving 



Mfc = Cu (ft • p,) iP, ■ Pk) / d\,,a 

Jo Jo 

^ g(ft • Pj,£)9iPk ■ Pj,e) 



dX 



{Xj^aXj.b) 



l-2e 



(23) 



The exchange diagrams contribute to two color struc- 
tures, symmetric and antisymmetric with respect to the 
generators Ta and on line Pj that are linked by gluon 
exchange to the lines Pi and Pk , respectively. Combining 
Fig. [TJd and c, we thus write 



Mb{Pi,e) + MaiPi, e) = M^^liPi,e) + M^llipi, e) 



(24) 



The first, antisymmetric, combination has the same color 
factor as the three gluon diagram. Fig. [T^, discussed 
above. It is given by 



5 



g{l3i ■ (3j,e)g{Pk ■ Pj,e) 



l-2e 



(25) 



where the difference in color factors in square brackets 
produces a commutator for the color matrices on the (3j 
line. Clearly, the integrals in parenthesis are identical 
and cancel. This result applies to arbitrary masses for 
the Wilson lines. As in the massless case, only the sym- 
metric contribution of the double exchange diagrams sur- 
vives and contributes to the soft matrix. As described in 
Ref. however, the symmetric contribution to the soft 
matrix at two loops is already generated by the exponen- 
tial of the one-loop anomalous dimension in Eq. ([3]). 

The reasoning above has a significant phenomenolog- 
ical application to 2 — > 2 production processes. Using 
Wick rotation, the vanishing of 3E diagrams after one- 
loop renormalization when pairs of invariants are equal 
applies to 2 — !■ 2 processes involving the production of 
pairs of heavy quarks from light quarks or gluons. In 
particular, we consider Wilson line velocities correspond- 
ing to momentum configurations with t ~ {pi — ps)^ — 
{Pi - PiY = M in Eq. ([J) with n = 2 16] . Note that 
because it is trivial to reintroduce /3|-dependence, the 
result applies as well to the limit where one or more line 
becomes massless. 

Of special interest are the anomalous dimension ma- 
trices that enter pair production: the 2x2 matrix for 
qq QQ and the 3x3 matrix for gg QQ. At 
one loop, these matrices are diagonal in their s-channel 
singlet-octet bases at m = t [9|. Specifically, the off- 
diagonal elements of these ■* are all proportional to 
ln{u/t). One- and two-loop diagrams contributing to 
these anomalous dimension matrices are illustrated by 
Fig. [3l where the single cikonal lines represent incom- 
ing light partons (quark pairs or gluons) and the double 
lines heavy quark pairs. At ii = i, the off-diagonal ze- 
ros in Fg directly reflect cancellation between pairs of 
diagrams such as those in Fig. [3^ and b. 

To illustrate the pattern of cancellation, we continue 
the notation above and write the amplitude, Mu corre- 
sponding to eikonal diagram D as the product of a color 
factor and a velocity-dependent factor, 

m},"' -Cz,4")(/5,). (26) 

In these terms, the vanishing of matrix elements that 
describe the mixing of color singlet-octet tensors in F^^' 
follows from identities relating the two diagrams, 

- F^l\M, (27) 

C2a = ~C2fc , (28) 




(a) (b) 




(c) (d) 




(e) (f) 

FIG. 3: One-loop 2E diagrams (a,b) that cancel ai u = t. (c) 
representative 3E diagram for this process. Diagrams (d,e) 
are two-loop diagrams that cancel. Diagram f represents the 
color structure of the webs. 



which are easily verified at one loop for u ^ t. We have 
shown above that at two loops with u = t, 3E diagrams 
like Fig.[3l: vanish independently of their color structure, 
leaving only 2E diagrams and the color-symmetric parts 
of 3E exchange diagrams at two loops. This suggests that 
at two loops the same pattern of cancellation may persist 

(2) . . . 

for T g . This is indeed the case, as we now explain. 

At one loop, the expansion of Eq. ([3]) in terms of one- 
loop anomalous dimensions can be thought of as the sum 
of contributions from each giuon-cxchange diagram. Ex- 
panding the one- loop form of Eq. ^ to order a^, we 
generate all color-symmetric contributions of two-loop 
ladder diagrams, corresponding for example to the color- 
symmetric contributions of Fig. [T]i and e above. The 
sum of all 2E diagrams is thus easily rewritten in terms 
of the expansion of Eq. ^ in terms of v'^g^ , plus two- loop 
contributions that involve commutators of generators on 
one or more of the Wilson lines. 

Examples of diagrams with antisymmetric color struc- 
ture, which do not correspond to the expansion of the 
one-loop anomalous dimensions, are Figs. [Tf and[3]i and 
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e, and the color-antisymmetric part of Fig.[T^. These dia- 
grams contribute only their antisymmetric combinations 

(2) 

of color generators to , and are thus proportional to 
Ca times a one-loop color factor. Figure[2]: is a schematic 
representation of the color structure of such diagrams and 
of their contributions to the two-loop anomalous dimen- 
sion matrices. 

We note that exactly as for the massless case in Ref. 
[m , all 2E diagrams are generated from an exponential of 
"webs" Webs are themselves 2E diagrams that are 
two-particle irreducible under cuts of the Wilson lines, 
starting with single-gluon exchange. The contributions 
of two-loop webs are precisely the antisymmetric color 
structures that we have just identified. 

As we have observed, Eq. holds for the color fac- 
tors of Fig. [3^ and b. The same relation then holds for 
the web contributions of Fig[3Ji and e, because they are 
proportional to the same color structure as the one-loop 
diagrams, as illustrated in Fig. [3]:. Therefore, whenever 
the identity for velocity factors, Eq. ([27]) . holds at two 
loops, this pair of two-loop diagrams cancels. The 2E 
diagrams of Fig. [3^ and b, however, can depend only on 
the invariants formed from the two eikonal velocities in 
question, so that at u = i, 

^ F^!\Pi) i^-t). (29) 

This relation will hold as well for diagrams associated 
with crossed ladders, as in Fig. [T}5, and other diagrams 
with color-antisymmetric contributions, including those 
with self-energies of the exchanged gluon and vertex cor- 
rections on the massive eikonal. As a result, cancellations 
between pairs of 2E diagrams that occur at one loop recur 
at two. For 2^2 kinematic configurations with u = t, 
then, we have in place of Eq. (jH) the relation 

rif(/3/)-i5(/3/)r^','(/3/) iu = t), (30) 

where D{j3i) is a matrix that is diagonal in the s-channel 
singlet-octet basis. For massless two-loop cases, the ma- 
trix D(f3i) is proportional to the identity, and we recover 
a special case of Eq. (U). More generally, however, the 
integrals of the 2E webs depend on masses. We will give 
explicit results elsewhere. In any case, the relation ([50]) 
applies for scattering at ninety degrees in the center of 
mass, including the limiting case of production at rest 
(s = 4mQ), of particular relevance to threshold resum- 



mation for the total cross section for heavy quark produc- 
tion . A subset of these diagrams have been analyzed 
very recently in [T9| . 

For massless particles, as noted in it is certainly 
natural to generalize the result Eq. (|4]) for lightlike par- 
tons to all orders. This possibility has been explored 
recently in Refs. [1^. The result that we have found 
here, that Eq. (|4]) does not apply to arbitrary massive 
kinematics, suggests that the zero-mass case is indeed 
special. We believe this should encourage the search for 
a symmetry or principle underlying the relation, Eq. 

In summary, we have shown that for products of mas- 
sive Wilson lines, one- and two-loop soft anomalous di- 
mensions are generally not proportional. We have noted, 
however, that diagrams that link two massless with one 
massive line through the three-gluon coupling vanish by a 
variant of the reasoning for three massless lines. We have 
also shown that diagrams that link a single massless line 
with two massive lines cancel for the case of two-to-two 
production processes at threshold, and when u = t, that 
is, for production at ninety degrees in the center of mass 
system. For these momentum configurations, the two- 
loop anomalous dimension matrix is diagonal in the same 
color basis as the one loop, although they are not related 
by a simple proportionality as in the massless case. The 
color-antisymmetric parts of exchange diagrams linking 
three eikonal lines cancel independently of masses. We 
carried out a numerical evaluation of the three-gluon dia- 
gram in Euclidean space, but the qualitative conclusions 
of this paper extend to Minkowski space. Whatever the 
field-theoretic origin of the results described above for 
massive and massless partons, the renormalization of the 
composite operators that link Wilson lines will be rel- 
evant to the analysis of light-parton jets, heavy quarks 
and potential new, strongly interacting particles at the 
Tevatron and the LHC. 
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